background = styleColorBar(iris$Petal.Length, ‘steelblue’), backgroundSize = ‘100% 90%’, backgroundRepeat = ‘no-repeat’, backgroundPosition = ‘center’
Nota: Esta es una aproximación cruda al modelo epidemiológico, y no debería tomarse como un resultado oficial ni de alta confiabilidad. Su interés es meramente referencial.
Usaremos la ecuación de Gompertz:
\[ y(t) = a \times e^{-e^{\frac{(\mu \times e)}{a} \times (\lambda - t) + 1}} \] Donde los parámetros a ajustar son a, mu (\(\mu\)) y lambda (\(\lambda\))
El resultado de usar todos los datos disponibles se ven a continuación.
Formula: y ~ a * exp(-exp(mu * exp(1)/a * (lambda - t) + 1))
Parameters:
Estimate Std. Error t value Pr(>|t|)
a 142770.74 262682.25 0.544 0.5918
mu 1360.23 2042.18 0.666 0.5117
lambda 47.44 22.24 2.134 0.0433 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 35.11 on 24 degrees of freedom
Number of iterations to convergence: 10
Achieved convergence tolerance: 7.901e-06
Además, esta ecuación tiene un valor de AIC de 273.6049097.
El gráfico que acompaña a este modelo, presenta (en rojo) la extrapolación a 7 días de la curva actual. Además, se muestra, en color naranja, la curva extrapolada usando los datos de los primeros 12 días desde el primer caso confirmado (que fue el 20202-03-06). Y como referencia, en en gris oscuro se muestra la curva estimada con los datos hasta las semana anterior.
Formula: y ~ a * exp(-exp(mu * exp(1)/a * (lambda - t) + 1))
Parameters:
Estimate Std. Error t value Pr(>|t|)
a 1104.6142 211.2713 5.228 6.81e-05 ***
mu 50.7101 2.8312 17.911 1.80e-12 ***
lambda 9.1077 0.4772 19.086 6.42e-13 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 16.74 on 17 degrees of freedom
Number of iterations to convergence: 119
Achieved convergence tolerance: 4.136e-06
Formula: y ~ a * exp(-exp(mu * exp(1)/a * (lambda - t) + 1))
Parameters:
Estimate Std. Error t value Pr(>|t|)
a 3935.43 5429.91 0.725 0.4870
mu 104.13 107.96 0.965 0.3600
lambda 14.66 6.70 2.188 0.0564 .
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 3.3 on 9 degrees of freedom
Number of iterations to convergence: 8
Achieved convergence tolerance: 5.963e-06
Esta visualización esta hecha usando el lenguaje de programación R, y hace uso de una serie de librerías (flexdashboard, ggplot2, echarts4r, DT), para producir los distintos gráficos y tablas.